Magneto-optical control of light collapse in bulk Kerr media 
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Magneto-optical crystals allow an efficient control of the birefringence of light via the Cotton- 
Mouton and Faraday effects. These effects enable a unique combination of adjustable linear and 
circular birefringence, which, in turn, can affect the propagation of light in nonlinear Kerr media. 
We show numerically that the combined birefringences can accelerate, delay, or arrest the nonlinear 
collapse of (2+l)D beams, and report an experimental observation of the acceleration of the onset 
of collapse in a bulk Yttrium Iron Garnet (YIG) crystal in an external magnetic field. 
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Waves propagating in multidimensional self-focusing 
media are subject to instabilities that lead to the catas- 
trophic collapse after a finite propagation distance, fol- 
lowed by beam filamentation or material damage [l|, H, 
i i, S, i 0, S i E El. In the (2+l)D [(2+1)- 
dimensional)] setting, above a certain threshold value 
of the input power, the critical collapse driven by fo- 
cusing nonlinear it ies is a universal scenario, observed in 
high-power excitations of plasmas Q, hydrodynami- 
cal systems [H , Bose-Einstein condensates (BECs) 0, Hj , 
and optical media [2I, 0, Q. In particular, in optical 
pulse propagation through amorphous media and crys- 
tals without special symmetries, the dominant nonlin- 
earity is the Kerr (cubic) effect, which is modeled by the 
nonlinear Schrodinger equation (NLSE) [l2[. Collapsing 
beams in Kerr media were studied in detail, especially in 
the course of the past decade [gI, 0, & & llO, El • 

The challenge to control the wave collapse (and in par- 
ticular to mitigate its detrimental effects) has recently 
drawn much attention [1, 0, [13, E| • As recently demon- 
strated, the collapse distance of ultra- intense laser pulses 
in air can be controlled by passive optical elements 
and in BECs the collapse time is strongly affected by the 
so-called Feshbach-resonance technique ^] . In condensed 
optical media, where the collapse occurs with pulse en- 
ergies far below the creation of plasma 0, 0] , a recently 
proposed scheme for collapse management relies on the 
use of a layered structure with the nonlinearity strength 
alternating in the longitudinal direction [10] (a similar 
mechanism was proposed for the stabilization of BECs 
in the 2D case [Uj). However, this scheme is difficult to 
implement in bulk media, and such structures may give 
rise to linear losses induced by the reflection of light from 
interfaces between the alternating layers. 

An alternative approach to control the transition to 
the collapse may be provided by optical birefringence, 
which promotes energy and phase transfer between the 
polarization components of the beam. The interplay be- 
tween birefringence and nonlinearity is known to induce 



coupling of the polarization rotation and a characteristic 
temporal evolution of solitons in optical fibers pj. In 
this paper, we explore the birefringence as a tool for the 
"management" of the collapse of (2+l)D beams, which 
may be relatively easily implemented in affordable ex- 
perimental conditions by using magneto-optical (MO) ef- 
fects jl^, El, El] 5 while avoiding reflection losses. To 
this end, we present numerical and experimental studies 
of the combined effects of linear and circular birefrin- 
gences on the dynamics of collapsing beams in a bulk 
self-focusing Kerr medium. We find that the onset of 
collapse can be accelerated, delayed, or even suppressed, 
at certain values of the combined birefringence strengths. 
We also show that the required linear and circular bire- 
fringence can be induced in a transparent MO Yttrium 
Iron Garnet (YIG, Y3AI5O12) crystal by the application 
of an external dc magnetic field, thus opening up the per- 
spective of using MO effects to generate and control var- 
ious nonlinear phenomena in optics. In our experiments, 
an adjustable balance of linear and circular birefringences 
was realized via a combination [14] of the Cotton-Mouton 
(or Voigt) and Faraday [16] MO effects in a bulk YIG 
crystal. Following the propagation of femtosecond pulses 
in the crystal, we observed a controllable decrease of the 
threshold input power necessary for the onset of collapse 
at the output facet as a function of the magnetic field. 
This also constitutes a first experimental study in non- 
linear optics where birefringence can be switched on and 
varied continuously in an adjustable fashion. 

The evolution of the complex electric-field amplitudes, 
Ur and ui^ representing the right- and left-circular polar- 
izations (RCP and LCP), in the presence of a Kerr non- 
linearity and combined linear and circular birefringences, 
obeys the coupled NLSEs, in the scaled form [l3| : 
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where z is the propagation axis, Vj^ is the transverse 
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Laplacian, while b and c are the strengths of the circular 
and linear birefringences, respectively. As evident from 
Eqs. ([1]) the linear and circular birefringences account 
for, respectively, the rates of linear amplitude mixing and 
phase shift between the RCP and LCP fields. In terms 
of the RCP and LCP, the ratio between the cross- and 
self-phase modulation coefficients (XPM/SPM) is 2 for 
cubically-symmetric crystals, including YIG [ll, [13]. 

Assuming vorticity-free solutions with circular symme- 
try, (2+l)D Eqs. ([1]) reduce to a (1+1)D form, in terms 
of z and the radial coordinate R. We consider input 
Gaussian beams, Ur{R,z = 0) = Aexp ^— cos^ and 

ui{R^z = 0) = Aexp (^—^^ sin6>, with normalized in- 
put width p=l, amplitude A, and symmetric polarization 
content 6>=7r/4. This choice corresponds to an unchirped 
Gaussian profile launched with an horizontal linear po- 
larization. Since the input does not carry vorticity, the 
circular-symmetric structure of the solutions is not sub- 
ject to an azimuthal modulational instability 0. 

Figure 1 shows propagation maps obtained by direct 
simulations of Eqs. ([1]). In the absence of birefringence 
[Figs. l(a)-(d)], the RCP and LCP components are cou- 
pled only via the XPM term, which becomes significant 
only close to the collapse point. While a low input power 
excitation [Figs. l(a,b)] results in beam diffraction, the 
collapse occurs after a finite propagation distance [for 
instance, z = 7.3 in Fig. 1(c)] when the input power ex- 
ceeds the critical level [8|] . A substantial converging por- 
tion of the phase fronts emerges near the collapse point, 
see Fig. 1(d). The collapse can be accelerated by the in- 
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FIG. 1: (Color online) Propagation maps of Ur{R, z) obtained 
from direct solutions of Eq. ([!]). (a)-(d): Zero birefringences, 
with (a),(b) low and (c),(d) high input powers. The pairs of 
panels (a),(c) and (b),(d) display the evolution of the inten- 
sity, \ur{R^z)\^ ^ and phase of Ur{R,z), respectively, (e), (f): 
Intensity maps in the presence of the birefringences. Numer- 
ical parameters are indicated above each respective panel. 




FIG. 2: (Color online) Evolution maps in the parameter 
space (6, c) (see the text), as obtained from direct simulations 
of Eqs. JI]), for A=1.13 and 0<z<50, up to ^ = 50 or the 
value Zcow < 50 in which the collapse occurs. The solutions 
that do not collapse up to z = 50 are represented in yellow 
areas. The dashed (green) curves show the calculated depen- 
dence between b and c corresponding to crystalline YIG in an 
external magnetic field, (a) ^coii, and (b) beam width at Zcoii- 



troduction of linear birefringence, as shown in Fig. 1(e). 
Circular birefringence, if acting alone, does not affect the 
collapse, since the respective terms in Eqs. ([1]) can be 
eliminated by a straightforward transformation. When 
both of the birefringences are present, the propagation 
distance necessary for the onset of collapse can be ex- 
tended^ i.e., the onset of the collapse is delayed [Fig. 1(f)] 
due to the interplay between amplitude and phase mix- 
ing. For low b and c values, such as those used in Fig. 
l(e,f) and typically achievable in experiment, the differ- 
ences in the evolution of Ur and ui are marginal, i.e., 
the RCP and LCP beam components feature the same 
collapse dynamics. With larger birefringence parameters 
the differences between the components become substan- 
tial; however, such large birefringence values were not 
accessible in the current experiment, see below. 

The results of systematic simulations are summarized 
in Fig. 2 by means of maps in the plane of (6, c). Panels 
(a) and (b) show, respectively, the values of 2:coii where 
the beam collapses, if Zcoii < 50, and the final values of 
the beam width, p{z = 50); in cases where the collapse 
occurs at Zcoii < 50, the final width is set as p = [black 
areas in (b)]. While Figs. 1 and 2 display the results for 
the RCP component, the LCP maps are similar in the 
entire range considered. As seen in Fig. 2, for a given 
input power the domination of amplitude mixing between 
RCP and LCP (c > b) usually results in an acceleration 
of the collapse, while dominant phase mixing {b > c) can 
lead to the delay or effective suppression of the collapse. 

In the experiment, a bulk YIG single-crystal was 
placed with the easy crystallographic axis [100] paral- 
lel to cf. Ref. These cubic dielectric crystals are 
highly transparent for optical signals in the near-infrared 
and exhibit large MO transmission coefficients, owing to 

1611. We chose to work at 



their ferrimagnetic phase lH, 
the wavelength of 1.2 /im, which offers an optimal trade- 
off between the magnitudes of the MO coefficients and 
the absorption losses [l^. The temporal dispersion in 
YIG is normal and weak in the near-infrared [l3], and 
hence a high power beam is free from the development of 
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FIG. 3: (Color online) Experimental setup schematics. OPA, 
optical parametric amplifier laser; T, cylindrical telescope; 
QWP, quarter wave plate; P, polarizer; ND, neutral density 
attenuators; LI, LI, aspheric coupling and imaging lenses; C, 
crystal sample. Insets: (a) Photograph of the magnet's inte- 
rior with lenses LI, L2 and crystal C; (b) The beam's waist 
profile at the input facet of the crystal. 



temporal modulational instabilitiess 

The application of an external magnetic field H per- 
pendicular to the propagation axis renders the crystal 
optically uniaxial, with the optical axis parallel to H 
and the respective linear birefringence proportional to H 
[3, [l5[. The largest phase retardation associated with 
the linear birefringence, for components of the wave vec- 
tor k perpendicular to the magnetic field, is 0.45 jam/ cm 
in YIG [lH at a saturation field of i7 = 800 G. This cor- 
responds to a saturation value of c = 0.1 in terms of Eq. 
(1) for our crystal. For wavevector components parallel 
to the magnetic field, the Faraday effect induces different 
refractive indices for the RCP and LCP waves (i.e., circu- 
lar birefringence), with a respective polarization phase re- 
tardation of 2 /im/cm at the saturation level [l6[. While 
the Faraday effect is predominant in the k || H geome- 
try it is weaker in the present setting, as only near 
the collapse point the wavevectors feature significant lat- 
eral components. An estimation corresponding to the 
case shown in Figs. 1(c), (d) yields an effective saturation 
value of b = 0.05. Below saturation, the Faraday coeffi- 
cient is proportional to the magnetic field, b ^ while 
the Cotton-Mouton coefficient depends on H quadrati- 
cally, i.e. c ~ 3, 15, 16|, implying a parabolic de- 
pendence between the birefringence coefficients in a given 
external magnetic field, b ^ c^. Calculated relations be- 
tween b and c for the YIG crystal at different magnetic 
fields are shown by the dashed curves in Fig. 2. 

The experimental setup is sketched in Fig. 3. We used 
a Spectra Physics model 800C Optical Parametric Am- 
plifier laser system, delivering pulses of 200 fs duration 
at a repetition rate of 1 KHz with peak powers < 50 
MW. The input beam was shaped by means of a cylin- 
drical telescope (T), followed by the combination of a 
quarter wave plate (QWP) and a polarizer (P), which 
were used to fix the input linear polarization parallel to 
H (i.e. horizontal). Variable neutral density filters (ND) 
were used to control the input power. The YIG sample 
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FIG. 4: (Color online) (a) Output beam profiles as a func- 
tion of the input peak power, in the absence of an external 
magnetic field, (b) Line-outs along x (full circles) and y (hol- 
low rectangles) around the central spot of the collapsed (2.9 
MW) output beam, (c)-(e) Output beam images, as function 
of the input power (horizontal) and the external magnetic 
field (vertical), for the fields: = 200 G (c), 400 G (d), 800 
G(e). 



was placed inside an electromagnet (GMW model 3470). 
The application of a driving voltage between the poles 
(±V) induced a uniform out-of-plane dc magnetic field 
in the crystal. A digital Tesla-meter (Group3 DTM-133) 
was used to calibrate the free-space magnetic field and 
verify its homogeneity in the sample. A combination of 
aspheric lenses, LI and L2, was used for coupling and 
imaging, respectively. Importantly, the lenses and crys- 
tal were mounted on nonmagnetic pyrex holders, see Fig. 
3(a). Standard metallic mounts, which usually hold mi- 
croscope objective lenses, were not used, as we found that 
they strongly affect the uniformity of the magnetic field 
within the sample. The sample input facet was placed at 
the focal plane of LI. The beam's profile at this plane is 
shown in Fig. 3(b). The FWHMs of this elliptical input 
beam were 30 and 15 /am along the major (x) and minor 
{y) axes. The beam profile at the sample output facet 
was imaged by L2 onto a Vidicon infrared camera. 

Experimental results are summarized in Fig. 4. Im- 
ages of the output beam as a function of the input peak 
power, in the absence of a magnetic field, are shown in 
Fig. 4(a). At the critical power of 2.9 MW, a central 
high-power spot appears in the diffracting background. 
This spot is circularly symmetric [see Fig. 4(b)], al- 
though the input beam profile was elliptical, with the 
above-mentioned major-to- minor axes ratio of 2 : 1. The 
spot fits a Townes profile [6], which indicates that it is 
generated by the onset of the collapse. The correspond- 
ing critical power complies with the theoretical predic- 
tion for collapse Q, using Aq = 1.2 //m, the linear re- 
fractive index of YIG no = 1.83 and its Kerr coefficient 
n2 = 7.2 X 10""*^^ cm^/W p^. The corresponding col- 
lapse distance prediction with A= 1.13, Zcoii = 7.3 [Figs. 
l(c),(d)], when rescaled back into physical units using 
the well-known transformations @] and assuming 

a circular input beam of 10 jam diameter, matches the 
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FIG. 5: (Color online) (a) Measured critical power required 
for the onset of the collapse in the output facet, as a func- 
tion of the applied magnetic field, (b) Numerically calculated 
squared critical amplitude of the input beam (A^r) necessary 
for the onset of the collapse at z = 7.3, as a function of the 
birefringences corresponding to crystalline YIG. (c) The cal- 
culated dependence between h and c in YIG, used in (b). 



sample length, 3 mm. At higher powers, filaments are 
observed in the output facet images, indicating that the 
collapse occurred earlier in the crystal. Following the ap- 
plication of a magnetic field, the onset of collapse at the 
output facet is observed at lower powers [Figs. 4(c-e)], 
with the largest collapse-acceleration effect recorded at 
a magnetic field of H = 400 G [see Fig. 4(d)l, corre- 
sponding to half the saturation field of YIG [17]. This 
observation agrees with the fact that, when YIG is ex- 
cited by light traveling perpendicular to the magnetic 
field, the linear birefringence is stronger than its circu- 
lar counterpart ESj] . A qualitative characterization of 
the output beam polarization state has also shown that 
without an external magnetic field the beam remained 
linearly polarized, while with an application of the field 
the beam became elliptically polarized, with the ellip- 
ticity growing with the field. This indicates the certain 
presence of magnetically-induced polarization dynamics 
in the crystal. 

Figure 5(a) shows the measured critical power, corre- 
sponding to the data of Fig. 4, that was required for the 
onset of collapse at the output facet, as a function of the 
magnetic field. In Fig. 5(b), this is compared to results 
obtained by the numerical solutions of Eqs. ([1]), where 
for each set of birefringence values (corresponding to a 
given magnetic field) we find the critical amplitude Acr 
of the input beam for which the collapse occurs at the 
propagation distance corresponding to the output facet 
{z = 7.3). Figure 5(c) again displays the relation between 
the normalized birefringence parameters of the YIG crys- 
tal used, cf. the dashed curves in Fig. 2. The behavior 
of the collapse detuning is similar in Figs. 5(a) and 5(b), 
even though the theoretical model did not take into con- 
sideration magnetically-induced losses. Specifically, the 
Cotton-Mouton effect is always accompanied by linear 
magnetic dichroism [15], and the Faraday effect entails 
circular magnetic dichroism [H, [l^, both of which are 
weak but present at A = 1.2 /im. 

In conclusion, we have investigated the combined ef- 
fects of circular and linear birefringences on the propa- 



gation of collapsing (2+l)D beams in self-focusing bulk 
Kerr media, and have shown that the onset of collapse 
can be accelerated, delayed, or suppressed, depending on 
the relative birefringence strengths. Experimentally, we 
have demonstrated a controlled acceleration of the col- 
lapse at the output facet of a ferrimagnetic YIG crys- 
tal, following the application of an external magnetic 
field which induces the birefringences. The direct ob- 
servation of magnetization-induced effects in collapsing 
beams provides a unique demonstration of an all-optical 
magnetically-controlled lensing mechanism, pioneering 
the use of MO crystals in nonlinear optics experiments. 
Finally, since Eqs. ([1]) are also Gross-Pit aevskii equations 
describing a binary BEG with linear interconversion [3] , 
accounted for by the coefficient c, similar phenomena may 
also be observed in nonlinear matter- wave dynamics. 
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